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ABSTRACT 

Aims. We performed a detailed pulse phase resolved spectroscopy of the accreting binaiy X-ray pulsar Her X-1 in the energy range 
3.5-75 keV and have established pulse phase profiles for all spectral parameters. For three parameters, the centroid energy of the 
cyclotron line, the photon index and the flux of the 6.4 keV iron line, we have studied the variation as a function of 35 d phase. 
Methods. We analyzed RXTE observations of the Main-On of November 2002 using data from the PCA and the HEXTE instruments. 
Four different time intervals of about ~1 d duration were selected to provide a good coverage of a complete Main-On. The intervals 
are centered at 35 d phases 0.03, 0.10, 0.15, and 0.20, respectively. 

Results. All spectral parameters show a strong modulation with pulse phase. While the centroid energy of the cyclotron line follows 
roughly the shape of the pulse profile, showing higher values close to the peak of the X-ray pulse, both the photon index and the iron 
line intensity exhibit distinct minima around the peak of the X-ray pulse. With respect to variations of the observed profiles with 35 d 
phase, we find that there is a clear evolution of the shape of the pulse profiles (flux versus pulse phase), a moderate increase of the 
maximum cyclotron line energy (found around pulse phase 0.7), but no significant evolution of the shape of the pulse phase profiles 
of the cyclotron line energy, the spectral power law index or the iron line intensity. 

Conclusions. The variation of spectral parameters as a function of the pulse phase provides important information about the system: i. 
the disappearance of the Fe line flux near the highest continuum flux may be an indication of a hollow cone geometry of the accretion 
structure; ii. The apparent non-dependence of the cyclotron line energy profiles on 35 d phase provides a new possibility to test the 
model of free precession of the neutron star, proposed to be responsible for the systematic variations in the pulse profiles. 

Key words, binaries: general - stars: neutron - X-rays: binaries - X-rays: individuals: Her X-1 - pulsars: individuals: Her X-1 



1. Introduction 

The X-ray source Her X- 1 was discovered in 1972 by Uhuru 
dTananbaum et al.lll972h . and classified as an accreting X-ray bi- 
nary. Her X-1 is one of the brightest and most studied persis- 
tent binary X-ray pulsars. The distance to the system is ~7kpc 
and the masses of the neutron star and its co mpanion are ap- 
proxi mately 1.5 Mq and 2.2 Mq, respectively (i Revnolds et al.l 
[1993). The X-ray flux shows periodic modulation on several 
different time-scales: pulsations with 1.24 s due to the spin of 
the neutron star, eclipses due to the 1 .7 d orbital period, and a 
35 d super-orbital period due to obscuration of the X-ray emit- 
ting region by the precessing accretion disk. This 35 d period- 
icity shows different states: two On-states with the ~11 d Main- 
On (at phase 0-0.31) and the ~5d Short-On (at phase 0.57 - 
0.79), separated by two Off-states (see e.g . Giacconi et apl 19731 
Gerend & Bovnton' '1976'; 'Bovnton et all ll 9801: IScott & Leahvl 



1999; .Scott et al.. .2000; Klochkov et alj l20()6h . The accretion 



disk is thought to be tilted (with respect to the orbital plane), 
warped and counter-precessing with a somewhat variable period 
of ~35 d. The onset of the flux (often identified with 35 d phase 
0.0) is called the Turn-On (TO) and corresponds to the transition 
from the Off-state to the Main-On. At this time, the outer rim of 
the disk clears the view to the X-ray emitting region close to the 
polar caps on the surface of the neutron star, leading to an in- 
crease in flux. Viceversa, the decline of the flux towards the end 
of the Main-On is identified with the inner edge of the accretion 
disk blocking our view to the X-ray emitting regions. A similar 



cycle is responsible for the Short-On with a maximum flux of 
roughly one third that of the Main-On. 

Her X-1 was also the first accreting X-ray pulsar for which 
a c yclotron line in the X-ray spectrum has been discovered 
dTriimper et al.ll9'78b . This absorption-like feature, now referred 
to as Cyclotron Resonance Scattering Feature (CRSF), is ob- 
served around 40keV and allows to estimate the neutron star's 
magnetic field. Applying the formula B 12 - (1+z) Ecyc/1 1.6keV 
(where B12 is the magnetic field strength in units of lO'^ Gauss, 
z is the gravitational redshift and Ecyc is the centroid energy of 
the cyclotron line), the first direct measurement of the magnetic 
field of a neutron star was achieved (~ 3 x lO'^ Gauss for Her X- 
1). The source shows a positive correlation between the cen- 
troid energy of the cyclotron line and th e bolometric luminosity 
(IStaubert et al.ll2007HVasco et al.ll201 Ih . 

Another observational feature of Her X-1 is the evolution 
of the pulse profile as function of energy and time, or better: 
35 d phase. For both the Main-On and the Short-On, a system- 
atic variation in the shape of the pulse profiles is found. These 
changes consist of the disappearing of some features in the pulse 
profile with an increase in energy and/or time. These variations 
repeat systematically on a time scale of ~3 5d. The se changes 
have been well known for some time (s ee e.g. lDeeteT et al. 19981 
Scott et al. 2000: Trum per et al.l ll986V Recently, 'sTaubeit et d] 
(2010a, b, 2012) on the basis of RXTE observations, gave a quan- 
titative description of these changes by providing a template, 
which can be used to predict the shape of the profile in the 9- 
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Table 1. Summary of the four intervals of the Main-On of 35 d 
cycle no. 323. The columns are: the number of the interval, the 
time interval in MJD, the exposure time, the corresponding 35 d 
phase interval, and the center of the 35 d phase interval. 



Interval 


limits 


exposure 


limits 


center 




MJD 


[ks] 


35 d phase 


35 d phase 


1 


52595.060 - 52596.743 


25 


0.007 - 0.056 


0.03 


2 


52597.959 - 52598.719 


25 


0.090 - 0.112 


0.10 


3 


52599.603 - 52600.726 


25 


0.137 - 0.161 


0.15 


4 


52601.316 -52602.406 


28 


0.186 - 0.217 


0.20 



13 keV energy range for any 35 d phase during the Main-On. An 
extension of this work for other energy ranges and including the 
Short-On is in preparation. 

[Triimper et al. (1986), on the basis of observations by 
Exosat, had proposed that free precession of the neutron star 
is responsible for those pulse profile changes, due to the vari- 
ation of the angle of the line of sight to the X-ray emitting re- 
gions on the surface of the neutron star. For Her X-1, the de- 
pendence of the spectral parameters on pulse pha se have been 
observ ed by severa l missions. Early ex amp les are Pravdo et a 1.1 
(11978 1) with OSO 8. lPravdo etaD(ll97 9') and'Soons et al.'('l990h 
with the A-4 experiment of HEAO-1 and McCrav et al. ( 1982) 
with HEAO 2 the Einstein Observatory, Voges et al. (|1982D 
with th e MPE/ AIT ballo n experiment, Kahabka (1987') with 
Exosat, Kunz et ap ( Il996h with Mir-HEXE. Later examples are 
lEndo et all (12000^ wi th ASCA , 'Zane & Ramsay" (2001) with 
XMM-Newton, Lutovi nov et alJ (2000) with the ART-P tele- 
scope on board of GRANADA and lKlochkov et alj (120081) with 
INTEGRAL. 

In this work we present a deep investigation of the variation 
of key spectral parameters as a function of pulse phase with the 
finest resolution so far and, for the first time, a discussion of a 
possible evolution with 35 d phase. Changes of the spectral pa- 
rameters with pulse phase are quite common among accreting X- 
ray pulsars and are generally attributed to a change in the view- 
ing angle of the accretion region (see e.g.. iKrevkenbohm et al.l 
[2004 and references therein). Here, we focus on the variation of 
the centroid energy Ecyc of the CRSF around 40keV, the pho- 
ton index T and the iron line characteristics. We also address 
the question whether our data support the free precession model 
which has been suggested to explain the variation in shape of the 
pulse profile. 



2. Observations and data analysis 

To carry out this analysis we used observations of Her X- 1 of the 
Rossi X-ray Timing Explorer (RXTE) performed in November 
2002. These observations correspond to cycle number 323 (cycle 
numbering according to Staubert et al. 2009b) centered at MJD 
-52599.32 with an observed Turn-On of MJD 52594.80. The 
data of this cycle provide the best coverage of any Main-On of 
Her X-1 observed with RXTE. Within this cycle, we sum up 
observations from four selected time intervals at different 35 d 
phases in order to have a good coverage of an entire Main-On, 
avoiding dips and eclipses. All intervals have a duration of about 
one day (details are given in Table [TJ. Figure [1] shows the four 
intervals selected for this analysis. 

For barycentric and binary corrections or bital param- 
eters and the ephemeris of IStaubert et al.l (l2009b were 
used: Porb= 1 -700167287 s, = +(2.8 ± 0.2) xlO"^^ ^ s"', a ■ 
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Fig. 1. PCA light curve of the Main-On of cycle 323. The grey 
zones represent the four intervals used for this analysis. For each 
interval the centered 35 d phase is given. Details of the time and 
35 d phase intervals are listed in Table [1] 
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Fig. 2. Mean 9-13keV pulse profiles of the four time intervals 
at diff'erent 35 d phases (see Table [TJ: Interval 1 in black (solid 
line). Interval 2 in green (dashed line). Interval 3 in red (dash 
dotted line) and Interval 4 in blue (dash dot dot d otted line). For 
general 35 d dependence see lStaubert et al.l (1201 Oaub) . As timing 
reference for pulse phase zero we use the "sharp-edge" at the 
trailing edge of the right shoulder of the main pulse which leads 
into an eclipse like minimum (Staubert et al. 2009b). Note that 
these pulse profiles are statistically so accurate that the uncer- 
tainties are smaller than the line width. 



sin /= 13.1831 It-s and 7^=52599.486440 MJD. To align the 
pulse profiles, a reference time to = MJD 52594.869900091 and 
a spin period f spi„ = 1.237761809 s were used. As "pulse phase 
zero" the sharp-edge feature at the trailing edge of the right 
shoulder of the main peak was used (see Staub ert et al.ll2009l) . 

The PCA was used in the energy range 3.5-60keV and 
HEXTE in the energy range 20-75 ke V. The data have been ana- 
lyzed with XSPEC 12. 6.(0 using the hi ghecut[3 spectral model. 



' http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec 

~ http://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XSmodelHighecut.html 
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Fig. 3. Cyclotron line centroid energy profiles as function of the pulse phase for Interval 1 (top left). Interval 2 (top right). Interval 
3 (bottom left) and Interval 4 (bottom right). 



power law continuum with an exponential cut-ofF. The differen- 
tial photon flux f(E) is described by: 



f(E)^A- 



, if£ <£cut 
£-r.exp(^), if£>£cut 



(1) 



where F is the photon index, Ecut is the cut-ofF energy and 
Efoid is the e-folding energy. In addition, two line features are 
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Fig. 4. The cyclotron line centroid energy as function of pulse 
phase for the four different 35 d phase intervals (see Table [T]i: 
Interval 1 in black (asterisks). Interval 2 in green (diamonds). 
Interval 3 in red (triangles) and Interval 4 in blue (squares). 



included in the model: a multiplicative Gaussian optical depth 
for the Cyclotron Resonant Scattering Feature (CRSF) and an 
additive Gaussian emission line for the iron fluorescence feature 
around 6.5 keV. T his model is the sam e one as used in earlie r 
analyses (see e.g.. ICoburn et al.ll2002l and IStaubert et al.ll2007h . 
We have consistently found that the highecut spectral model 
describes the spectrum of Her X- 1 best. The discontinuity of this 
model at the cut-ofF energy does not constitute a problem and we 
have verified that no systematic uncertainties are introduced in 
estimating the spectral parameters. We used the PC A data up 
to 60 keV, which has become feasible by new response matrices 
(PCA response vll.7, 2009 May 11). These high energy data 
can definitively improve the photon statistics at energies around 
and beyond the cyclotron line featur e. The same choice h as been 
used in previous analyses ( see e.g . Ferrigno et al.ll20lTI) and it 
was confirmed by Rothschil d et al.l ( 201 Ih with RXTE observa- 
tions of Cen A. The recorr0 spectral component was also added 
to the model to normalize the background. After we had ver- 
ified, that the energy and width of the iron K line is constant 
(6.45 keV and 0.5 keV, respctively), we fixed those values in or- 
der to minimize the number of free parameters. GoodXenorfl ob- 
servational modes of the PCA for high resolution analysis were 
used to exctract light curves and spectra. An additional system- 
atic uncertainty of 1% was added in using the PCA data. To an- 
alyze the data, version v6.11 of the HEASOFT software^ and 



http://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/ 
XSmodelRecom.html 

http://heasarc.gsfc.nasa.gov/docs/xte/abc/ 
pca_issues.html#configs_modes 
^ http://heasarc.gsfc.nasa.gov/ftools/ 
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Fig. 5. The cyclotron line profiles of Fig. 3 for the four 35 d intervals (black crosses): Interval 1 (top left). Interval 2 (top right). Interval 3 (bottom 
left) and Interval 4 (bottom right), normalized to a common amplitude. The big red filled dots represent the weighted mean values of all four 
intervals. The solid blue line represents a best fit function (a combination of two cosine components). We conclude that, within the given statistics, 
there is no evidence for any change in shape of the profiles with 35 d phase. 



background estimation files Sky_VLE dated 2005 November 28 
were used. We have investigated whether there are correlations 
between the fit parameters: no dependence of any single param- 
eter on the others has been found. All spectral results stated be- 
low have been obtained through spectral fits with the assumption 
of negligible absorption due to neutral hydrogen, both intrin- 
sic to th e source and interstellar (which is at 5.1 x 10'^ cm"^, 
[Dal Fiu me et al.l[T9 98). The introduction of absorption as a free 
parameter does not lead to a reduction in x^- 

For the pulse phase resolved analysis, we have extracted 
the spectra in 14 pulse phase bins over the 1.24 s pulse, the 
widths of which were chosen to have roughly identical photon- 
statistics in each spectrum. The centers of the 14 phase bins 
are at pulse phases: 0.1, 0.3, 0.45, 0.55, 0.625, 0.675, 0.7125, 
0.7375, 0.7625, 0.7875, 0.8250, 0.8750, 0.9250, and 0.9750. For 
the cyclotron line energy E^yc this is the finest resolution in pulse 
phase ever achieved (the four smallest bins around the peak of 
the pulse have a width of 1/80 of a phase). 

With the spectral function described above, good fits with re- 
duced between 0.9 and 1 .2 are generally achieved. However, 
splitting up the data into small pulse phase intervals reveals that 
Isager x^ values appear: up to ~12 for the pulse phase range 
0.775-0.85, and 1.4-2.8 for the pulse phase range 0.5-0.775. 
The main reason for this is the existence of the so called "10 keV 
feature" (see the discussion in 3.5). Modeling this feature by an 
extra Gaussian is generally successful and brings the x^ down 
to acceptable values (in a few cases is still found around 1.3, 
which is most likely due to an imperfect modeling by a simple 
Gaussian). 



3. Pulse phase variability 

We concentrate on the detailed study of the variation of the fol- 
lowing spectral parameters as function of pulse phase: the cen- 
troid energy Ecyc of the cyclotron line, the photon index F, and 
the intensity of the iron line at 6.4 keV. The phase profiles of 
these spectral parameters are compared with the shape of the flux 
profiles (the "pulse profiles") for the four different 35 d phase in- 
tervals. In Figure |2] the 9-13keV pulse profiles of the four in- 
tervals are shown. Note that the pulse profiles are of such high 
statistical quality that the uncertainties of the flux values (in each 
of the 128 pulse phase bins) are comparable to the width of the 
lines chosen to plot the profiles. The profiles are repeated in the 
following figures (as dashed lines) to allow a direct comparison 
with the profiles of the spectral parameters. For the remaining 
spectral parameters, such as the cut-off energy Ecut, the fold- 
ing energy Efoid (with reference to the used spectral model - 
highecut), and the width and depth of the cyclotron line (tTcyc 
and Tcyc), we will restrict ourselves to profiles of the combined 
data of cycle 323. For these parameters the statistics is such, that 
no statements can be made about a dependence on 35 d phase. In 
all Figures, the uncertainties given are one sigma (68%) values. 

3.1. Cyclotron line energy 

The profiles of the cyclotron line energy E^yc as a function 
of pulse phase are shown in Figures |3] and |4] Generally, Ecyc 
roughly follows the shape of the pulse profile: the broad maxi- 
mum is found to be close to the peak of the pulse profile (around 
pulse phase 0.7) for all 35 d phases. The formal uncertainties 
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Fig. 6. Photon-index profiles as function of pulse phase for Interval l(top Jeft), Interval 2 (top right), Interval 3 (bottom left) and Interval 4 
(bottom right). 



of Ecyc (as determined by XSPEC) are of the order of a few to 
several percent for pulse phases 0.0-0.6 and 0.9-1.0 and around 
one percent or less for pulse phases 0.6-0.9 (close to the peak 
of the pulse). Adopting an additional systematic uncertainty of 
O.lSkeV (~ 0.4%) for all data points, the profiles are well rep- 
resented by cosine functions (see below). We find slightly dif- 
ferent values for the mean Ecyc and the peak-to-peak ampli- 
tude: a slight, but significant increase is found as function of 
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Fig. 7. The photon index F as function of pulse phase for the 
four different 35 d phase intervals (see Table [T}: Interval 1 in 
black (asterisks). Interval 2 in green (diamonds). Interval 3 in 
red (triangles) and Interval 4 in blue (squares). 



35 d phase. The weighted mean values of the four highest val- 
ues of Ecyc of each group are 40.16 + 0.22 keV, 40.89 + 0.20 keV, 
41. 02+0. 15 keV and 41. 52+0.20keV for mean 35 d phases 0.03, 
0.10, 0.15 and 0.20 (intervals 1 through 4), respectively. This 
corresponds to an increase by ~ 0.7 ke V per 0. 1 units in 35 d 
phase. The minimum Ecyc values have larger uncertainties with 
a tendency to decrease with 35 d phase, but, within uncertain- 
ties, they are consistent with a constant value (around 33.3 keV). 
The peak-to-peak amplitude is then slightly increasing with 35 d 
phase (with a mean of ~ 7.6 keV). 

While the values for the maximum Ecyc and the peak-to-peak 
amplitude are dependent on 35 d phase, the sinusoidal shape of 
the modulation does not change. In order to verify this statement, 
the following analysis was performed: the original Ecyc profiles 
(Figure HJi were scaled such that the separation of minimum to 
maximum Ecyc had a constant value of 7.64 keV, using the min- 
imum values measured in phase bin 0.2-0.4 and the weighted 
mean of the four highest values in each of the four profiles as 
reference. In this way, normalized Ecyc profiles were generated. 
These normalized profiles are shown in Fig. |5] each in compar- 
ison to the weighted mean values and a smooth best fit function 
of the weighted mean: a combination of two cosine functions 
with common mean (37.31 keV), amplitude (3.82 keV), period 
and phase-zero values of 0.70 and 0.772 for pulse phases 0.3- 
0.875 and 1.47 and 0.3 for the remaining phases, respectively. 
Visually, there is no obvious deviation of any of the four data 
sets from the weighted mean. As quantitative statistical test we 
performed the following: comparing the individual normalized 
profiles with that profile representing the weighted mean leads 
to values (normalized to the degrees of freedom) of 0.72, 
1.46, 0.72 and 1.32 for the four intervals, respectively. As an 
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Fig. 8. Iron line intensity profiles as function of pulse phase for Interval l(top left), Interval 2 (top right). Interval 3 (bottom left) and Interval 4 
(bottom right). All intensities less than 1.0 are consistent with zero, we therefore show two sigma upperlimits. 



additional test, we have then performed a Monte Carlo simula- 
tion, in order to answer the following question. Assuming the 
profiles are indeed the same for all four intervals, that is equal 
to the best fit mean profile: what is the probability to find an in- 
dividual profile which deviates from the mean profile in such a 
way that the reaches or exceeds the observed value? We have 
simulated 100 profiles for each of the four intervals, reproduc- 
ing the statistical conditions of each interval (that is, the width 
of the Gaussian distribution from which the random numbers 
were drawn were chosen according to the uncertainty of the cor- 
responding measured data point). Those probabilities are 68%, 
29%, 73% and 4% for the four intervals, respectively. Even 4% 
for Interval 4 is a rather high and acceptable probability. We con- 
clude, that, within our statistical accuracy, we find no evidence 
for any change in the shape of the E^yc profiles with 35 d phase. 

3.2. Photon-index 

Figure |6] shows the profiles of the photon index F as function of 
the pulse phase. This spectral parameter is relatively constant for 
pulse phases up to 0.6 (showing a small positive slope), but has 
a distinct dip around the pulse peak and a smaller one around the 
left shoulder. The photon-indeces range between 0.42 + 0.01 and 
1.07 + 0.01 for Interval 1, between 0.54 + 0.01 and 1.08 + 0.02 
for Interval 2, between 0.36 ± 0.01 and 1 . 10 ± 0.49 for Interval 3 
and between 0.33 + 0.01 and 1.07 + 0.01 for Interval 4. In gen- 
eral, the spectrum gets harder when the flux increases. While the 
main features of the pulse phase dependence of F are very much 
the same in the four 35 d intervals, the depth of the new narrow 
dip around phase 0.62 shows a clear trend towards smaller values 
from Interval 1 to Interval 4 (see Fig.|7). 



3.3. Iron line intensity 

Another spectral parameter which shows a clear modulation as 
function of pulse phase is the iron line intensity (see Figure [8]l. 
It shows a moderate increase until the left shoulder of the pulse 
(phase ~ 0.7), after which it drops almost to zero around pulse 
phase 0.8, before reaching the start value again at pulse phase 
0.9. The minimum is coincident with the peak of the X-ray pulse 
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Fig. 9. The iron line intensity profiles for the four different 35 d 
phase intervals (see Table [T]i: Interval 1 in black (asterisks). 
Interval 2 in green (diamonds). Interval 3 in red (triangles) and 
Interval 4 in blue (squares). 
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Fig. 10. Pulse phase profiles for the cutoff energy Ticut i^op ieft), the folding energy Efoi^ {top right), the width of the cyclotron line c {bottom 
left), and the depth t^^^ of the cyclotron line {bottom right) for the total data sets. 



(see Figure [8] and |9|l in a similar way as the photon index. The 
iron Hne intensity ranges between (for all intervals) and 12.40 
for Interval 1, 14.62 for Interval 2, 13.19 for Interval 3 and 10.75 
for Interval 4. The profiles of 35 d intervals 2^ are quite simi- 
lar to one another, while for Interval 1, the iron line intensity is 
significantly lower for pulse phases < 0.6. 

3.4. Other spectral parameters 

For the other spectral parameters Ecut, Efoid the width cr^yc and 
the depth Tcyc, we show their pulse phase dependence only for 
the total data sets (not resolved in 35 d phase), because in these 
parameters we find a considerably larger scatter, such that no 
conclusions can be drawn about a possible variation (or not) 
with 35 d phase (the 35 d phase resolved profiles can be found 
in IVasco 20120. The profiles are shown in Figures [10] Also 
for these parameters there is significant modulation with pulse 
phase, and they generally show a broad maximum around the 
phase of the pulse peak, with the interesting feature that for Efou 
there is a dip in this maximum, right at the pulse peak. It is also 
worth noting, that, except for Efoid, the measured values for pulse 
phases < 0.4 are rather high, not following the shape of the pulse 
profiles. 

3.5. The "10 keV feature" 

It has been noticed (e.g. ICobum et al.ll2002h . that in fitting ob- 
served spectra of accreting X-ray binaries with a cut-off power 
law a "wiggle" or "bump" around ~ 10 keV to ~ 15 keV is found 
in the residuals. This applies to many different sources includ- 
ing Her X-1 (also to non-cyclotron line sources), to observations 
with different satellites (e.g. Ginga, Beppo/SAX and RXTE), and 



* http://tobias-lib.uni-tuebingen.de/frontdoor. php?source_opus=6346 



to analysis with different spectral functions involving a power 
law (e.g., highecut ,plcut, Fermi-Dirac or NPEX). The fea- 
ture is therefore considered to be inherent to accreting X-ray bi- 
naries. So far, no physical explanation has been suggested. 

In our current analysis, using data from observations of 
Her X-1 by RXTE, this feature does not play a role in the phase 
averaged analysis. However, in the pulse phase resolved analy- 
sis, larger values appear due to the existence of this "10 keV 
feature": up to ~12 for the pulse phase range 0.775-0.85, and 

1.4- 2.8 for the pulse phase range 0.5-0.775. Modeling this fea- 
ture by an extra Gaussian is generally successful and brings the 

down to acceptable values (in a few cases is still found 
around 1.3, which is most likely due to an imperfect modeling 
by a simple Gaussian). The centroid energy is generally found at 
~ 16 keV with a cr of ~ 5 keV. Fig. [TT]shows an example of the 
spectrum of Interval 3, pulse phase bin 0.80-0.85, where the x^ 
is reduced from 5.55 to 1.38 for 224 dof. We have verified, that 
all the other spectral parameters are not altered when introduc- 
ing this extra Gaussian. We note that this feature occurs most 
strongly at a pulse phase close to the peak of the pulse, where 
all other spectral parameters also tend to show their strongest 
variability. 

4. Summary and discussion 

We have selected the best data available in the RXTE archive 
of an observation of a Main-On of Her X-1 to perform de- 
tailed pulse phase resolved spectroscopy in the energy range 

3.5- 75 keV. The observations of November 2002 ( cycle no. 323 ) 
provide the best RXTE coverage of a Main-On (IVascoll2012l) . 
such that an attempt could be made to study pulse phase depen- 
dent spectral variations also as a function of the phase of the 35 d 
modulations - flux and shape of the pulse profiles -, which are 
believed to be due to precessional motion of the accretion disk 
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and, possibly, the neutron star itself. As spectral model we have 
used a power law with an exponential cut-off (the highecut 
model in XSPEC) modified by two line features, a multiplica- 
tive cyclotron line (in absorption) and an additive iron fluores- 
cence line (in emission). We have concentrated on three spectral 
parameters: the centroid energy of the cyclotron line, the pho- 
ton index T of the power law component (valid between 3 keV 
and ~ 20keV), and the intensity of the iron fluorescence line at 
6.4 keV. For these three parameters, the analysis was performed 
for four one day long intervals at different mean 35 d phases 
(0.03, 0.10, 0.15 and 0.20, respectively). For the other spectral 
parameters of the applied spectral model, no 35 day phase de- 
pendence was studied. 




5 10 20 50 

Energy (keV) 



Fig. 11. Residuals of the X-ray spectrum of pulse phase bin 
0.80 - 0.85 for Interval 3 without (top) and with (bottom) an extra 
Gaussian component in the fit. 



4.1. Pulse phase dependence 

The results of the pulse phase resolved analysis can be summa- 
rized as follows. We distinguish between results which confirm 
and improve in statistical significance which was known before 
and new results. 

4.1.1. Confirmed and improved results 

- All spectral parameters show a strong modulation with pulse 
phase. 

- The strongest variations always appear around the peak of 
the pulse (at maximum flux). 

- The profile of the centroid energy Ecyc of the cyclotron line 
follows roughly the modulation of the pulse profile (maxi- 
mum energy around maximum flux). The shape is close to 
"sinusoidal". 

- The photon index profile (Fig. |6]and|7]l shows a strong dip 
close to the pulse peak. The dip around the main peak is a 
known feature, nicely confirmed here and very well resolved 
in this analysis. 

4.1 .2. New results 

- It is the first time that for Her X-1 a phase profile of the 
centroid cyclotron line energy Ecyc is produced with such 
a high resolution and statistical accuracy. The four small- 
est bins around the peak of the pulse have a width of 1/80 
of a phase. The mean variation (minimum to maximum) is 
~ 23% (Figs. [3] and |4|. The modulation is close to "sinu- 
soidal"; however, at least two sinusoidal components (with 
greatly different periods) are needed to model the modula- 
tion. 



- For 35 d intervals 1-3, the peak of the Ecyc profile is struc- 
tured (see Fig. |3]l, with local maxima near the maximum of 
the main peak of the pulse profile and around the position of 
the left shoulder of the main peak. A similar structure is seen 

in Tcyc. 

- The folding energy Efoid follows the pulse profile, except for 
a narrow dip in the peak of the pulse (Fig. [10] upper-right). 

- Similarly, the cut-off energy Ecut, the width and depth of the 
cyclotron line (iXcyc and Tcyc), roughly follow the pulse pro- 
file, except for higher values at phases 0.0-0.5 (Fig. [TOt . 

- While the strong dip in the photon index profile (Fig. |6] and 
|7]) close to the main pulse peak was known (even though less 
well measured), the second, very narrow dip (around phase 
0.61) is a newly detected feature. We note that this is close to 
the position of the left shoulder component of the main peak. 

- The profile of the intensity of the 6.4 keV iron line show a 
distinct minimum at maximum flux (Figs. [8] and|9]l. 

- We find that the so called "10 keV feature" seen in the 
spectra of several accreting binary pulsars, is around 16keV 
in Her X-1, with a strong pulse phase dependence: essen- 
tially concentrated in the narrow pulse phase interval 0.8- 
0.9. Spectral fits for all the other phase intervals do not show 
significant residuals. 

4.2. Precession phase dependence 

With respect to any variation of the spectral parameters (or their 
pulse phase variations) with 35 d phase, we consider only the cy- 
clotron line energy Ecyc, the photon index F and the intensity of 
the iron fluorescence line at 6.4 ke V. For the other parameters the 
photon statistics is not sufficient to distinguish between different 
35 d phases. It is the first time that a study of the variation of 
spectral parameters with 35 d phase was done. So, the following 
results are all new. 

- The maximum value for Ecyc (reached near the pulse peak) 
shows a slight, but significant increase with 35 d phase. The 
mean value is 41 .0 + 0. 1 keV, and the slope is ~ 0.7 keV per 
0. 1 units of 35 d phase. 

- The minimum value for Ecyc (around pulse phase 0.3) is for- 
mally consistent with a constant value of 33.5 + 0.5 keV (we 
would not rule out a small decrease). Together with the max- 
imum increasing, this means that also the peak-to-peak am- 
plitude (mean value 7.6 + 0.5 keV) is increasing with 35 d 
phase. 

- The shape of the Ecyc profile (see Figs. |4] and |5]l is indepen- 
dent of 35 d phase. 

- Similarly, the shape of the F profiles is independent of 35 d 
phase, except for the new narrow dip around pulse phase 0.62 
(Figs.|6]and|2l), where a trend to a smaller depth from Interval 
1 to Interval 4 is seen. 

- The shape of the iron line intensity profiles are similar to one 
another, except for Interval 1 which shows a distinctly lower 
intensity for pulse phases <0.6 (Figs.[8]and|9]l. 

4.3. Comparison with previous results and physical 
conclusions 

For Her X-1, the dependence of the spectral parameters on pulse 
phase have been studied by several authors on the basis of ob- 
servations with different X-ray instruments (see Introduction). 
In the following we highlight our new findings and significant 
improvements and discuss possible physical implications. 
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Pulse phase dependence of the centroid cyclotron line en- 
ergy: This modulation is quite common among accreting X-ray 
pulsars and is generally believed to be due to the changing view- 
ing angle under which the X-ray emitting regions are seen (e.g., 
iKrevkenbohm et al] |2004| and references therein). If we adopt 
the idea, that peaks (high X-ray flux) in the pulse profile can be 
associated with beamed radiation emitted from the accreting re- 
gions around the magnetic poles, then the observed pulse profiles 
are a representation of the emission characteristics. During a full 
neutron star rotation we may be seeing different emitting spots 
(most likely from a complex magnetic field configuration which 
deviates from a simple dipole (see e.g. lSuchv et al.l l2008). possi- 
bly having a multi-pole configuration) and the emerging beamed 
radiation under changing angles. The situation is further compli- 
cated by gravitational bending. It is therefore not straightforward 
to interpret observed pulse profiles in ter ms of accretion geom - 
etry and beaming characteristics (see e.g. lCaballero et al"1l201 ll) . 
If we, however, assume that the main peak in the pulse profile of 
Her X- 1 is due to a reasonably narrow pencil beam (following 
e.g., iPravdo et alll 19781: iKlochkov et alj|2008h . emitted perpen- 
dicular to the neutron star surface at the hot spot of the main 
accreting pole, then we may associate this pulse phase with a 
situation where the observer is looking down into the accretion 
column, parallel to the magnetic field lines and the velocity of 
the in-falling material. At this phase we see the maximum cy- 
clotron line energy and the hardest spectrum (a sharp minimum 
in photon index F) - that is, we are looking deep down to small 
distances from the neutron star surface where the magnetic field 
is strongest and the temperature is highest. Here, we also find 
the largest width and depth of the cyclotron line, as well as the 
highest Ecut and a local dip within a broad peak of Efoid. 

Along similar lines, the newly discovered structure of a 
second (very sharp) minimum in F around pulse phase 0.61 
(Fig. |7ll, corresponding to the left shoulder of the main peak, 
may be the signature of a second pole. At this pulse phase we 
also find indications for another relative maximum in Ecyc (Fig. 

m. 

Pulse phase dependence of the 6.4keV iron line flux: Figs. 
[8] and |9] show a sharp and deep minimum of the intensity of the 
6.4 keV iron line, again coincident with the main pulse peak. The 
profile is very similar to that of F (Fig.|7ll and the Fe line intensity 
drops to near zero at pulse phase ~ 0.79 (where F is minimal). 
It is the first time, that such a sharp mini mum in the 6.4 keV 
flux is observed. While 'Choi et al.l d 19941) . in observations by 
Ginga, and lZane"e t al. ( 2004), in observations by XMM-Newton, 
had found a sinusoidal-lik e modulation w ith a broad minimum 
around the pulse peak. Oosterbroek et al.lfcOOO ). with data from 
Beppo/SAX, were only able to measure the Fe-L line around 
1 keV, for which they had found a minimum around the pulse 
peak which was similarly sharp and deep as we now find for the 
6.4 keV line. It is generally believed, that the origin for the Fe 
lines is fluorescence in relatively cool material which is being il- 
luminated by the X-ray beam. For the geometrical site, however, 
where this fluorescence is taking place, or whether it is predomi- 
nantly in reflection or in transmission, no general consensus has 
been reac hed so far. Severa l different regions have been proposed 
(see. e.g. IChoi et al.l [19941 and references therein): 1) the outer 
edge of the accretion disk extended in the vertical direction; 2) a 
corona of the accretion disk; 3) the atmosphere of the compan- 
ion star; 4) a surface at the Alfven radius. The Alfven surface is 
currently the most popul arone. Howe ver , also the accretion col - 
umn was proposed (e.g.. lLeahvll200ll and linam & Bavkalll2005h . 



despite the argument by I White etaP (Il983h that the material in 
the column may be too hot. 

We like to stress the observed fact, that at the pulse phase 
where we see the highest continuum flux there is essentially no 
Fe line flux. This should provide some constraints in finding a 
feasible geometry. If the above adopted scenario is correct, by 
which we look down into the pencil beam at the pulse peak, the 
Fe line result may be consistent with a hollow cone geometry, 
where the majority of the radiation escapes through the empty 
center of the cone, where there is no material which can flu- 
oresce, while the walls are seen at off-center angles. It is, of 
course, assumed that the bottom of the cone close to the neu- 
tron star surface is filled with X-ray emitting plasma. 

The non-dependence of the cyclotron line energy proflle on 
35 d phase and the model of neutron star free precession: One 
of the results of this study is the finding that the shape of the 
cyclotron line energy profile (E^yc vs pulse phase) is indepen- 
dent of 35 d phase (Figs.[3]and|4]i, while the maximum cyclotron 
line energy increases slightly by ~0.7 keV per 0. 1 units of 35 d 
phase. If we attribute the strong variations in pulse profile with 
35 d phase to precession of the neutron star, that is to the change 
of the angle under which we see the X-ray emitting regions, we 
would expect also the shape of the cyclotron line energy profile 
to change with 35 d phase. However, this expectation needs a 
quantitative analysis. Recently concluded detailed model calcu- 
lations (Postnov et al. 2012), assuming neutron star free preces- 
sion, have shown, that under a specific set of assumptions (in- 
cluding a certain multi-pole geometry and emission into narrow 
pencil beams), the observed pulse profiles in the e n ergy range 
9-13 keV (the same as used in'Staubertetal ]|20T0allbll201 2) can 
be reprod uced with high accur acy. We point out here that the 
model of iPostnov et alj (1201 2|) is an attempt to give a physi- 
cal explanation of the variation in pulse profiles by neutron star 
free precession, requiri ng several spec i fic a ssumptions, while the 
template approach of Staubert et al.l ('2012), mentioned in the 
Introduction, does not need any physical assumptions, it just 
makes use of the observed pulse profiles. We plan to investi- 
gate whether the here described constancy of the shape of the 
cyclotron energy profile is consistent with the physical model. 
In this sense our result opens a new channel to test the model of 
free precession of the neutron star 
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